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Abstract 

We analyze the optimal unambiguous discrimination of two arbitrary 
mixed quantum states. We show that the optimal measurement is unique 
and we present this optimal measurement for the case where the rank of 
the density operator of one of the states is at most 2 ("solution in 4 di- 
mensions"). The solution is illustrated by some examples. The optimality 
conditions proved by Eldar et al. [Phys. Rev. A 69, 062318 (2004)] are 
simplified to an operational form. As an application we present optimality 
conditions for the measurement, when only one of the two states is de- 
tected. The current status of optimal unambiguous state discrimination 
is summarized via a general strategy. 
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1 Introduction 

Among the subtleties in quantum information processing and in quantum com- 
munication protocols are the properties that originate from the fact that in 
quantum mechanics non-orthogonal states cannot be discriminated perfectly. 
In the most naive approach to quantum state discrimination - the minimum 
error discrimination (cf Ref [HE]) - this leads to the fact, that the identifica- 
tion of a state might be erroneous with some finite probability. Ivanovic [3] and 
Dieks [U showed that one can avoid erroneous measurement results and that a 
measurement with a conclusive state identification is possible. In the case of 
non-orthogonal states, this strategy cannot work with a success probability of 
one. Peres showed in Ref. [5] how the optimum of this success probability can 
be achieved in the case of pure states, both having the same a priori probability. 
The discussion of the optimal unambiguous discrimination of two pure states 
was completed by Jaeger and Shimony in Ref. [B]. They derived the optimal 
solution for arbitrary a priori probabilities. 

Although it was long ago stated to be an interesting problem [7] , the unam- 
biguous discrimination of mixed states did not attract much attention for a long 
time. This changed with an example introduced by Sun et al. in Ref. [S] and 
the first general analysis of the unambiguous discrimination of mixed states by 
Rudolph et al. in Ref. 0. After that, several general results and special classes 
of optimal solutions were found, cf. Ref. El [12 [131 HH US HI]. While 
Bergou et al. derived in Ref. [10) I12j the optimal measurement for the unam- 
biguous discrimination of a pure state and an arbitrary mixed state, no analysis 
so far did succeed to produce a general solution for the simplest instance of gen- 
uine mixed state discrimination, the discrimination of two mixed states where 
both density operators have a rank of 2. Also the simple question whether the 
optimal measurement in general is unique remained unanswered. 
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The answers to these two questions are among the central results of this 
contribution. The uniqueness of the optimal measurement is stated in Propo- 
sition [Tl] and the general solution for rank 2 density operators is presented in 
Sec. El 

A valuable tool to approach both questions turned out to be a result by 
Eldar et al. in Ref. [17j . They showed necessary and sufhcient conditions for 
a given measurement to be optimal. However, these conditions are difficult to 
verify, since the criterion implies the proof of the existence or non-existence of 
an operator with certain properties. In Corollary [5] we reformulate this criterion 
in such a way, that it can be directly applied to a given measurement. As a 
further immediate consequence of this Corollary we will be able to provide simple 
optimality conditions for a very special type of measurement: The measurement 
which only detects one out of the two states, cf. Sec. 15.11 It will also become 
possible to provide a simple proof and a deeper insight into the fidelity form 
measurement [IHIIII], cf. Sec. 15.21 

Before we arrive at these results, we first provide an analysis of unambiguous 
state discrimination (USD), beginning in Sec. [21 where we derive general results 
and continuing in Sec. [31 in which we specialize to the optimal case. 

An analysis of the structure of the optimal measurement in particular yields 
Theorem [H This Theorem is a cornerstone in order to prove the uniqueness 
of the optimal measurement and also provides a simple proof of the "second 
reduction" shown by Raynal et al. in Ref. [TT]. We summarize and deepen the 
analysis carried out in Ref. [TT] in Proposition [31 Proposition [SI and Lemma [3 

In Sec. [7] we will provide a generic scheme in order to approach a given 
optimization problem for USD. We conclude in Sec. [HI 

2 Defining properties of USD 
2.1 Main definitions 

In quantum state discrimination of n quantum states it is usually assumed that 
the density operators pi, . . . , p„ of all possible input states are known, together 
with the probability pi,...,p„ of their occurrence. For 1 < /i < n, the a 
priori probability > and the corresponding density operator > with 
tr(pp) — 1 naturally combine to a weighted density operator — PfiPfj.- Hence 
the trace of a weighted density operator 7^ is the a priori probability of the state, 
tr(7^) = p^. Using this notation, the input states are represented by a family 
of positive semi-definite operators S — (7^). For a meaningful interpretation 
in terms of probability, we clearly need to have X]/jtr(7^) = 1. However, we 
will not require this normalization, as the subsequent definition and analysis is 
independent of it, and for certain statements (cf. e.g. Proposition [2]) it will be 
useful to explicitly allow tr(7Ai) < 1- 

In the following we will only consider the case of two input states, i.e., 
/i = 1,2. We restrict our analysis to finite-dimensional quantum systems, such 
that any possible quantum state of the system can be represented by a density 
operator which acts on a Hilbert space of finite dimension. We will use 
the formalism of generalized measurements in which a physical measurement 
with M possible outcomes is described by a positive operator valued measure 
£ — {El, . . . , Em) on Jf, i.e., by a family of M positive semi-definite operators 
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which sum up to the identity, Ek ~ ^■ 

Let us introduce our notation. We denote by ker^ = € | = 0} 
the kernel of an operator A, and we write AJ^ — {A\(f>) \ |(/)) S Jif} for its 
image. The support of a positive semi-definite operator p is written as suppp = 
e \ 3a > 0: p - a\cj)){(j)\ > 0}. Note, that the support of p is the 
orthocomplement of its kernel, suppp = (kerp)"*- and since p is self-adjoint, 
pJ^ = supp p holds. 

By a projector we always mean an orthogonal projector, unless we explicitly 
state that the projector is oblique (cf. Lemma fTHl in Appendix E]) . We use upper 
case Greek letters for orthogonal projectors, = S = S^. The symbols "c" 
and "D" are used such that they also include equality, i.e., £/ = ^ if and only 
if ^ C and £/ Z) ^. 

For a pair of weighted density operators S — (71, 72), we abbreviate 

supp 5 EE supp (71 + 72 ) = supp 71 + supp 72 , ( 1 ) 

for the collective support of S, which is the physically relevant subspace for the 
discrimination task and keriS for the common kernel of S, which then is the 
trivial subspace, 

ker5 = kcr(7i + 72) = ker7i H ker72. (2) 

The task of optimal unambiguous discrimination of two mixed states is defined 
as follows. 

Definition 1. A positive operator valued measure £ = {Ei, E2, E-/) is called 
an unambiguous state discrimination (USD) measurement of a pair of weighted 
density operators S — (71,72) */tr(i?27i) = and tr(i?i72) = 0. The success 
probability Psucc of £ of S is given by 

-fsucc {£;S)=tr{Eiji)+tr{E2j2)- (3) 

A USD measurement £ of S is optimal if it has maximal success probability, 
i.e., if for any USD measurement £' of S, Psucc{£',S) > Psucc{£'',S) holds. A 
USD measurement £ of S is called proper i/supp(£'i + E2) C supp 5. 

The condition tr(_E27i) = is equivalent to suppi?2 C ker7i and tr(i?i72) = 
is equivalent to supp Ei C ker 72 . Thus it is simple to write down some USD 
measurement for a given pair S. In the next section we will see, that it is 
sufficient to consider proper USD measurements. But the set of proper USD 
measurements in particular is compact (this follows from the above definition 
or more directly from Proposition [3|) and hence there always exists at least one 
proper USD measurement, which maximizes the success probability. 

2.2 Trivial subspaces 

For any USD measurement £ = {Ei, E2, E?) of 5 = (71,72) one readily con- 
structs a proper USD measurement £' — (E[, E2, E!,) with the same marginal 
probabilities, i.e., tr(i?i7i) = tT{E'iji) and tr(i?272) = ti'(£;272)- For that the 
most straightforward approach is to choose E[ and E'2 to be the projection of 
El and E2 onto supp 5 and to set E!, — 1 — E'^ — E!^. 
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As an important feature of proper USD measurements we will show that 
the optimal proper USD measurement is unique (cf . Proposition llip . Such a 
statement of uniqueness clearly can only hold if we require that the measurement 
is proper. For illustrative reasons let us provide an example of an optimal USD 
measurement, which is not proper and where the measurement operators do not 
even commute with the projector onto suppiS: We consider two non-orthogonal 
pure states with 

7i-^|l>(l|, 72 = i|+>(+|, (4) 

where |+) = (|0) + |l))/\/2 and the measurement £ = {Ei,E2, 1 - Ei - E2) 
with 

i?^ = (3-3/V2)|e^)(e^|, (5) 

where |ei) = (|0) - |1) - |2))/V3 and jea) = (%/2|0) + |2))/^/3. It is straightfor- 
ward to verify, that this measurement is a USD measurement and has a success 
probability of Psucc = 1^1/ V2 as given by the optimal solution due to Peres 

The subspace ker5 cannot play any role in USD, since the support of 71 
and 72 is orthogonal to this space. Similarly, the subspace supp7i n supp72 
necessarily is orthogonal to the support of Ei and E2 , since supp Ei C ker 72 and 
supp E2 C ker 71 . The following proposition is a consequence of this observation: 

Proposition 2 (cf. Theorem 1 in Ref. lllj). Let S = (71,72) be a pair of 
weighted density operators. Denote by Iljj. the projector onto (ker 71 -|- ker 72) 
and write — (nj|7inj|., Xlji 72X1^1) for the projected pair. Let q>0. 

Then £ is a proper USD measurement for S with Pgncci^', S) ~ q if and only 
if £ is a proper USD measurement for with Psucci£',S^) — 9- 

(Note, that (ker 71 -I- ker 72) is the orthocomplement of (supp 7x0 supp 72), which 
can be considered to be the "parallel" part of the support of 71 and 72.) 

Proof. If f is a USD measurement of S, we have n||-£'^n||- — E^ and so clearly 
tr(i?^7i,) = tr(£^^nj|7i,n||) holds. We have that supp i?^ C supp S and supp E^ C 
(ker 71 -I- ker 72). Due to supp = supp S n (ker 71 + ker 72 ) , it follows that £ 
is also proper for 5*. 

For the converse, since £ is proper for 5* we have in particular IIifE'^IIij. — 
Ef^ and hence tr(£'^n^7yn^) = tr(i?^7,y). Furthermore we have suppE"^ C 
suppiS'''^ = supp 5 n (ker 71 4- ker 72), i.e., £ is proper for S. □ 

2.3 The role of E? 

For the discussion of USD measurements it is useful to note that the measure- 
ment operator corresponding to the inconclusive result, E-?, already completely 
determines a proper USD measurement. 

Proposition 3. For an operator Ei and a pair of weighted density operators 
S — (71,72) there exist operators Ei and E2, such that £ — (Ei, E2, E->) is 
a proper USD measurement of S , if and only if E-r acts as identity on ker 5, 
E? >0, 1- E-> >0 and 71(1 - E9)-f2 = 0. 

Given E-i, the proper USD measurement £ of S is unique. 
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Proof. It is straightforward to see that the conditions are necessary. The proof 
of sufheiency and uniqueness is constructive: Let us write Qi for the obhque 
projector from ker72 fl suppiS to supp7i n (ker7i + ker72) (for a brief intro- 
duction to obhque projectors cf. Lemma fTHl in Appendix!^. Then we have for 
any proper USD measurement E2Q1 = and Ei = EiQi. Hence 

El = Q\{Ei + E2)Qi = Q1(1 - E,)Qi (6) 

is the only candidate for Ei, given E?. Due to 1 — E? > 0, this construction 
ensures that i?i > 0. An analogous construction holds for i?2. 

It remains to show that Ei+ E2 — — E7) =0. We decompose the Hilbert 
space into the sum 

= ker5 © (supp7i + supp72). (7) 

With IIj^ the projector onto ker5, we have £'^11^ = and since E-/ acts as 
identity on ker5, also (1 — E->)Il± = holds. Using, that by construction 
7ii?^72 = 0, we furthermore have 

ji[Ei + £^2 - (a - £?)]72 = -71(1 - £^?)72 = 0. (8) 

From 7i(l — £7)72 = and 1 — £? > it follows that with IIn the projector 
onto supp7i n supp72, we have (1 — £?)n|| — 0. Furthermore, one verifies that 
Qiji = (1 - n||)7i and hence 71 (£1 + £2)71 = 7i£i7i = 7i(l - £?)7i- A 
similar argument for 72 finishes the proof. □ 

Due to this Proposition [3] we sometimes refer to an operator E9 as a proper 
USD measurement if it satisfies the conditions of the Proposition. In addition 
from the Proposition it follows easily that the set of USD measurements is 
bounded and closed and hence in particular it is compact. 

A proper USD measurement is already uniquely defined by £'7(72 — ji)E-? 
(as it will turn out below, cf. Lemma [TUl in the optimal case this operator is in 
some sense much simpler than £? itself). Namely, with n_L the projector onto 
keriS and (71 +72)" denoting the inverse of (71 +72) on its support we have 
the identity 

E? = IIx + (71 + 72)" {7172 + 7271 

+ Vtia/ V7i[72 - £7(72 - 7i)£^?]\/7i VTi 

^, (9) 

+ %/72 Y %/72[7i - ■E^?(7i - 72)£?]V72\/72 

}(7i +72)"- 

In order to see this, first note that using -^71(1 ~ £7)^/72 = and £7 > the 
term in curly brackets can be rewritten as 




Then due to (71 + 72) (71 + 72) = 1 — n_L and once more 71 (1 — £7)72 = we 
see that the right hand side of Eq. ® is given by 

n± + (71 +72)^(71 +72) [1 - (1 -£?)](7i +72)(7i +72)" 

= n_L + (i-ni)£?(i-ni) (n) 



6 



This expression is equal to E-?, since for a proper measurement E-rll^ = Il± 
holds. 

Using the forthcoming Lemma [TUl Eq. and Proposition [31 it will become 
possible to reconstruct the optimal measurement given only the projective part 
of E7. This projective part is given by ker(l — E?). It has a very specific 
structure, which originates in the condition 7i(l — i?7)72 = 0. Let denote the 
projector onto supp7^ and n_L denote the projector onto kerS. For any proper 
measurement these projectors satisfy ni(l — i??)n2 = and (1 — E?)!!^ = 0, 
and hence Lemma [T71 (Appendix E]) applies, i.e., for any proper measurement, 

ker(l - E->) = ({ker(l - E?) supp7i} 

+ {ker(a-£'7)nsupp72}) ®ker5 (12) 

holds. (Note, that in the right hand side of Eq. (|12p the first and second term are 
in general not orthogonal and share supp7i n supp72 as a common subspace.) 
Although this result may seem to be quite technical, in certain situations it 
turns out to be a quite powerful tool. 

3 Simple properties of optimal measurements 

The following theorem makes a simple but fundamental statement about the 
structure of optimal measurements. It states that apart from trivial cases no 
vector which is in the kernel of 71 or in the kernel of 72 will be in the support 
of E? . This clearly gives an upper bound on the rank of E-/ . On the other hand 
the condition 71(1 — E?)j2 = provides a lower bound on the rank of E-?. The 
second part of the theorem states that these bounds coincide and fix the rank 

of £;?. 

Theorem 4. Let £ = {Ei, E2, Ef) be an optimal USD measurement for a pair 
of weighted density operators S ~ (71,72)- Then (supp _£? Hker 71) = (supp iJvH 
ker72). 

If £ in addition is proper, then suppi?? n ker7i — ker5, suppi?? H ker72 = 
ker5, and ranki?? rank 7172 + dimker5. 

(Remember, that the rank of an operator A is given by dim{Ajf) = dim Jf^ — 
dimker^, i.e., the number of strictly positive eigenvalues of A^yl.) 

Proof. Let € suppS? n ker7i. Then due to G suppS? there exists an 
a > such that E-f — a\(j)){(j)\ > 0. We define a new USD measurement by 
£' = (£'i,i?2 + a\(j)){(j)\, E? — a\(j)){(j)\) . From the optimality condition for 8, 
i.e., Psucc{£',S) < Psucc{£,S), we find a(0|72|0) < which only can hold if 
72I0) = 0. Since S suppi??, (suppi?? n ker7i) C (suppi?? n ker72) follows. 
An analogous argument holds for the "d" part and finishes the proof of the first 
assertion. 

From this result by intersection with (kcr 71 ) one immediately finds (supp E-rD 
ker7i) = (suppi?? H ker5). In the case of a proper measurement, however, 
suppi?? D ker5 and hence (suppi?? n ker7i) = ker5 follows. 

Let E-, denote E-? projected onto supp 5. Since the measurement is proper, 
E? — E!, is the projector onto ker5 and suppi?? = suppi?,' © ker5. From the 
previous results we have i?7^nker72 = {0} and i?.'72^nker7i = {0}. Then 
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due to LcmmafTOl f Appendix if follows ker(72i?7) — ker El, and ker(7ii?7 72) = 
ker(iJ^72). Hence, 

dimkerii'.' = diniker(72i??) = dim ker(i?772) ^^^^ 
= diniker(7ii5772) — diniker(7i72), 

where we used that dimker A = dimker A^ for any operator A and that 71(1 — 
-©7)72 -71(1- £^?)72 = 0. □ 

3.1 Orthogonal subspaces 

An important consequence of the first part of Theorem [5] is the following 

Lemma 5. Let £ — [Ei, E2, E-?) be an optimal USD measurement for a pair 
of weighted density operators S — (71,72)- Suppose that 11 zs a projector with 
C (ker7i n supp5). 
Then Eill — if and only if E2TI = 11. 

Proof The "if" part follows directly from < UE->n = -UEiU. For the 
converse we have suppi?? D E^IiJ^ = (11 — E2U)Jf C ker7i and thus due to 
Theorem 21 E-fllJff C ker S. But since ker 5 is orthogonal to HJ^, we have 
IIE9II = 0. Thus = E-rU = H - E2TI. □ 

In particular let S2 denote the projector onto ker 71 n supp72. Then neces- 
sarily for any USD measurement i?i E2 = and hence by virtue of Lemma [3 
for any optimal measurement -E2S2 — S2 holds. With Ei denoting the pro- 
jector onto ker 72 H supp7i we obtain Ei'Ei ~ Ei in an analogous way. These 
observations are at the core of the following 

Proposition 6 (cf. Theorem 2 in Ref. |TI]). Let S — (71,72) be a pair of 
weighted density operators. Denote by Ilgkew the projector onto (ker7i-t-supp72)n 
(ker 72 -I- supp7i) and write 5""^°"" = (nskcw7inskow, nskew72nskow) for the pro- 
jected pair. Let E? and Ef^'^'" be two operators satisfying £^7*^°"" = E-? + {TL — 

Il-skcw ) . 

Then E-? is an optimal and proper USD measurement for S, if and only if 
Ef^'^^ is an optimal and proper USD measurement for S^^™ . 

In this case, the failure probability for £ of S is the same as for £^^™ of 

^skcw 

tr(7i + 72) - Psucc{£,S) = tr[n,kcw(7i + 72)] - Psucc(f 5^*^™). (14) 

(Note, that (ker 71 + SUPP72) H (ker 72 -I- supp7i) is the orthocomplement of 
(Si + S2) . For the projected pair S^^^"" , the spaces supp(nskow7inskow) and 
supp(nskow72nskow) are skew, where two spaces £/ and are called skew, if 

£/n^^ = {0} = i^n^-L.) 

Proof. Due to the discussion leading to the Proposition, for any optimal mea- 
surement E-/ we have i^vHskew = E-> and hence nskow(l — -E7)nskow = Ti--Ef^™. 
It follows that 1 - Ef'''^ > and that nskcw7inskcw(l - ^^-f ™)nskcw72nskcw = 
0. Furthermore we find due to Hgkcw^ ^ ker 5 that 

kgr^skcw ^ ^ 72)nskcw] = kern,kcw © ker 5, (15) 
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where both terms in the direct sum are orthogonal. This shows that Ef^°™ 
acts as a projector onto ker5^'*''". Since obviously Ef^™ > we have shown 
that Ef^™ is a proper USD measurement for S^^™. The converse, namely that 
E-? ~ Ef^'^'" — (1 — Ilskcw) is a proper USD measurement of S, in fact holds 
for any proper measurement Ef^^"" of S^^°™. This follows from Eq. (fT5|) and by 
noticing that for S^^"^ = {Ef'=^, £:f ™, ™) the measurement defined by E? 
is given by £ = (Sf ™ + Si, ™ + S2, £:?). 

In order to show that given E->, the measurement Ef^""" is optimal, suppose, 
that Ef^™' is proper and has a higher success probability than Ef^™. Then 
it is easy to see that E!, — Ef^^""' — (1 — Ilskow) would yield a higher success 
probability for S than i??, in contradiction to the assumption. 

On the other hand, since iJ? Ilskcw = E-?, any optimal and proper E-? mini- 
mizes tr(£'7nskcw(7i + 72)nskow)- But this is minimal for optimal Ef^'^'" , since 

77' t^skcWTT I — I 

Proposition [2] and Proposition [H] can be used independently from each other, 
in contrast to the original result in Ref . [TT] . Proposition [2] and Proposition El 
provide a method to obtain all optimal measurement^^ for a given pair S by 
considering a different pair S' where dimsupp5 > dimsupp5'. This is in 
particular useful, if dimsuppiS' < 4, since in Sec. [5] we will provide an analytical 
solution for any such pair. If dimsupp5' < 2, then the general solution can 
already be obtained due to the result by Jaeger and Shimony [B] . Also the pair 
S' might possess a two-dimensional common block diagonal structure which 
was not present in the original pair S and allows a solution of the problem (cf. 
Ref. [IH]; for a simple criterion in order to detect such structures, cf. Ref. T^). 
Apart from that, using both propositions all optimal measurements can be found 
by just considering pairs of states which do not possess any orthogonal (like 
supp7i nker72) or parallel (supp7i n SUPP72) components. 

The following property simplifies actual calculations. 

Lemma 7. With the notations of Proposition\^and Proposition\^let denote 
the (non-linear) mapping from S to and analogously Tskew the mapping from 
S to 5"'^™. 

Then Tjt ° T]t = T(|-,. ^skcw ° Tskow = Tskew and Tskow ° T-jf = r||. O Tskcw- 

Proof. We abbreviate t[5]^ for 7^ after the application of t, i.e., (t[5]i , t[5]2) = 
t[S]. One verifies 

supp Tj/f[S]f, = supp 7^ n supp n^i, (16) 

and 

supp Tskcw [5]^ = supp(l - 1:^)7^ (1 - S^) ^^^^ 
= supp 7^ n kerE^. 

From the first equation we immediately get kerT||[iS]i -|-kerr|f[5]2 = J^, i.e., 
Tj/f is acts as identity on T|j^[5]. In order to show that Tskcw is idcmpotent one 
verifies that suppTskcw[5]i n ker Tskcw['5]2 = {0}. 

Due to suppT|f[5]i n kerT||[5]2 = supp7i n ker72 it follows that 

(Tskcw O T|f)[5] = (S71S, S72S), (18) 



^ In the original work it was only shown that one may choose the measurements in 
that specific way. Here we showed that all optimal measurements must have this structure. 
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where^ = nskcwllij. = Hifllskow Analogously due to ker Tskow [5] i+kerTskow [5] 2 = 
ker 71 + kcr 72 we have 

O Tskcw)[5] = (S71S, S72S), (19) 

and thus the thhd assertion holds. □ 

As an important consequence one can apply the mappings r in any order 
and in particular due to (rgkcw ° — ''"skew ° Tit? ^ second application of both 
mappings is never necessary. 

The action of r[| on 5 = (71,72) is non-trivial, if and only if rank(7i +72) < 
rank7i + rank 72. Similarly, the action of Tskow is non-trivial if and only if 
rank 71 > rank 71 72 or rank 72 > rank 71 72. We call a pair of states S strictly 
skew, if (Tskow o ^)[S] = S. 

Let us briefly mention a convenient way to construct the mapping Tskcw ° t^- 
As shown in the proof of Lemma [71 we can write 

S' = (r,kcw o ^)[S] - (S71S, S72S), (20) 

with S = 1 — — El — S2. Now let (|sii)) and (|s2j)) be Jordan bases (cf. 
Appendix IC]) of supp7i and supp72. Then IIn = J2iex ^^"^ ~ 

EfcsD;,, \Sf,k){s^k\, with X = {k \ {sik\s2k) = 1}, = {« | Vj: {su\s2j) 0}, 
and3^2 = {j I Vi: (si,|s2j) = 0}. 

Summarizing Proposition [2] and Proposition [HI if £' = {E[,E'2,Ei,) is an 
optimal and proper USD measurement of S' = (rgkcw ° ■^n-)['5], then £ = {E[ + 
Til, E'2 + 1^2, E'^ — Y^i — S2) is an optimal and proper USD measurement of S. 
The optimal success probability computes to 

Psucc(f;5) =P,ucc(f';5')+tr[(Ei + E2)(7i+72)]. (21) 
3.2 Classification of USD measurements 

We want to introduce a classification of the different types of optimal measure- 
ments for USD. Given the dimension of supp5, the classification is according to 
the rank of the measurement operators. For a Hilbert space of dimension d, we 
consider the optimal and proper USD measurements £ = {Ei, E2, E?) for pairs 
of weighted density operators S = (71, 72). We restrict the analysis to the case, 
where Tgkow and r||- act as identity on S, i.e., to the case of strictly skew pairs. 
Then rank 7172 = rank 72 = rank 71 = r and dim ker 5 = d — 2r holds. All 
optimal measurements with ranki^i = ei and rank_E2 = £2 will be considered 
as one type of measurement, denoted by (61,62). As we will see in subsequent 
sections, the construction method of the optimal measurement mainly depends 
on the type of the measurement. The symmetry of USD for exchanging the 
label of 71 and 72 makes it only necessary to develop a construction procedure 
for the case where e.g. 61 < 62. Thus a measurement class [a,b] with a < b 
denotes both measurement types (a, b) and (6, a). We now count the number of 
measurement types and measurement classes. 

Since we consider proper measurements, we have supp Ei n supp E2 = {0} 
and hence 61 -f 62 = rank(i?i +E2) and < r. Let us denote by S the dimension 
of the projective part of E?, i.e, S = dimker(]l — E-/). Then ei + 62 + S — d and 
5 < ranki??. From Theorem |4] we have that rank£^? = ?■ + (d — 2r). On the 
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other hand, at least ker(l — E-?) D keriS, i.e., (5 > — 2r. In summary we arrive 
at the eonstraints 

ei < f , 62 < r, and r < ei + 62. (22) 

From the situation where 71 and 72 have a two-dimensional block diagonal struc- 
ture, one can see that for any possible ei and 62 which satisfy the constraints 
in Eq. (|22p . one can find a pair S = (71, 72) such that an optimal measurement 
is of the type (ei, 62). 

Counting the possible combinations to satisfy the conditions in Eq. (|22p . one 
finds 

#types= i(r + l)(r + 2) and #classes = [(§ + 1)^J , (23) 

where ^types denotes the number of measurement types and ^classes the num- 
ber of measurement classes. Here we used the floor function, \_x\ — max{fc G 
Z I fc < a;}. 

Measurements of the type (ei , 62) with 61+62 — r actually are von-Neumann 
measurements. (Obviously there are always r + 1 such measurement types.) 
This can be seen, since then d — r = rankiJ? >5 — d — ei — e2 = d — r, i.e., 
ranki?? = dimker(l — E-?) and hence Ej is projective. But then trE'i +tri?2 — 
tr(l — i??) = r = €1 + 62 holds and hence all eigenvalues of Ei and E2 are either 
1 or 0. This proves the assertion. 

As we will see in Sec. 15. H and Sec. 15.21 an analytic expression for the optimal 
measurement is only known for the class [r, r] and the special von-Neumann class 
[0,r]. These classes may occur for any r > 1 and thus in particular solve the 
two-dimensional case (r = 1) and "half" of the four-dimensional case (r — 2). 
The remaining two classes (one of which is von-Neumann) in four dimensions 
are solved in Sec. 16.11 and Sec. 16.21 

4 The optimality conditions by Eldar, Stojnic &; 
Hassibi 

Eldar, Stojnic, and Hassibi provided in Ref. [171 necessary and sufficient condi- 
tions for the optimality of a USD measuremeno 

Theorem 8 (Eldar, Stojnic & Hassibi [E]). Let £ ~ {Ei, E2, E-?) be a proper 
USD measurement for a pair of weighted density operators S — (71,72)- Denote 
by Ai the projector onto ker72 H suppiS and by A2 the projector onto ker7i H 
suppiS. This measurement is optimal, if and only if one can find an operator Z 
such that for 11 — 1,2, 

Z > 0, ZE? = 0, (24a) 
A^(Z-7^)A^ > 0, and A^(Z - 7^)^^ = 0. (24b) 

In Ref. [13, this statement was only proved for the case ker5 = {0}. How- 
ever, the generalization presented in Theorem [8] follows immediately from the 
original statement. 

^Indeed Eldar ei al. proved conditions for the optimality of a USD measurement for an 
arbitrary number of states. 
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In Theorem [S] necessary and sufficient conditions for optimality where pre- 
sented. However they are not operational, as the existence or non-existence of 
Z is difficult to prove. We show in Appendix |B1 that the unknown operator Z 
can be eliminated, and the above conditions can be re-expressed as follows: 

Corollary 9. With the 'preliminaries and notations as in Theorem\^ a proper 
measurement £ of S is optimal if and only if 

(Ai - A2)-B7(72 - 7i)£^?(Ai As) > (25a) 
(Ai - A2)E.{j2 - 7i)^^?(]l - -Bv) = 0. (25b) 

The conditions for an optimal USD measurement are now expressed as a 
series of equations and positivity conditions on only i??. Remember the fact 
that E-i already completely determines a USD measurement (cf. Proposition [3]). 

The first condition in the above Corollary [U Eq. ([25a|, relies on the fact, 
that a positive semi-definite operator in particular has to be self-adjoint. Thus, 
the condition in Eq. (j25ap is only a compact notation for the three conditions 

Ai£;?(72 - 7i)-E7Ai > 0, (26a) 
A2£;?(7i -72)£^?A2 > 0, (26b) 
Ai£;?(72 - 7i)£;?A2 = 0. (26c) 

(Obviously these conditions are sufficient for Eq. ()25ap . The necessity follows 
from multiplication of Eq. ()25ap by from the left and Qy from the right. 
Here are the oblique projectors as defined in the proof of Proposition [31 ) 

The second equation, Eq. (|25b|) . in Corollary [5] makes a statement about the 
projective part of Ei. This is the content of the following 

Lemma 10. Let E-? be an optimal and proper USD measurement of a pair of 
weighted density operators S = (71,72) with supp7i n supp72 — {0}. Denote 
by H? the projector onto supp-E? and by A the projector onto ker(l — £'?). 
Then £'7(72 — 7i)-E? = n7(72 — 7i)H7 = A(72 — 7i)A. 

Proof. Let H_l denote the projector onto ker5. Then due to supp 71 n supp 72 = 
{0}, we have (Ai - A2)Jf = (1 - H^)jr. Due to Ii^E7-/f, = 0, the optimality 
condition in Eq. (|25bp hence reads -£7(72 — 71 )£'?(! — £^7) = or 

£^?(72 - 7i)^^? = E->{j2 - li)E7^. (27) 

For the first equality we multiply this equation from the right by the inverse (on 
its support) of E-? and in a second step from the left and obtain the equations 

E?{-f2 - 7i)n7 = E->X-i2 - 1i)Et. (28a) 
nv(72 - 7i)H7 = n7(72 - 7i)i;7. (28b) 

Since the right hand side of the first equation is self-adjoint, the assertion follows. 

For the second equality we have E-r/S. = /S. and N = iJ?(l — A) = £'? — A with 
iVA = 0. Thus E-;{\-E-7) = N{\-N). But ker(l-A^) = ker(a-i;7+A) = {0} 
and hence the optimality condition in Eq. (|25bp reads Ei\^2 — ~ 0. Thus 



£^7(72 - -/i)E7 = £;7(72 - 7i)A (29a) 
A(72 - 71 )^^? = A(72 - 71) A (29b) 

holds, where in the second step we multiplied the first equation by A from the 
left. □ 



12 



Lemma [TU] is the key to prove the uniqueness of the optimal and proper 
USD measurement, since due to the identity in Eq. Q we have seen that any 
USD measurement is solely defined by -E?(72 — 7i)-E?- Hence in the case of 
supp7i n SUPP72 = {0}, the optimal and proper USD measurement can be 
uniquely determined, given II?, the projector onto the support of E-?. But 
since the set of optimal and proper USD measurements of S is by virtue of 
Proposition [5] equal to the set of optimal and proper USD measurements of <S*, 
having supp7i O supp72 — {0} (cf. also Lemma [7]), it remains to show that 
the support of E-/ is unique. This follows from the fact that the rank of E-/ is 
fixed by virtue of TheoremlH together with the convexity of optimal and proper 
measurements. Namely, for any two optimal and proper USD measurements E-? 
and E-?, also ^{E-> +£'?) is an optimal and proper USD measurement. But since 
E-f and E-/ are positive semi-definite, rank(£'7 + E-/) — ranki?? = rankiJ? can 
only hold if supp E? — supp E-f . Thus we have proved the following 

Proposition 11. For a given pair of weighted density operators, there exists 
exactly one optimal and proper USD measurement. 

5 Two special classes of optimal measurements 
5.1 Single state detection 

For certain pairs of weighted density operators S — (71, 72) it may be advanta- 
geous to choose a measurement with tr(_Bi7i) = 0, e.g. if tr(7i) is much smaller 
than tr(72). We refer to this situation as single state detection of 72. In the clas- 
sification scheme proposed in Sec. 13. 2( the single state detection measurements 
can be identified with the class [0,r] (where r = rank 7172). 

For a proper measurement, tr(i?i7i) = can only hold if already Ei = 0. 
If the measurement is optimal then due to Lemma [Sj Ei ~ implies E2 = A2. 
(The projectors were defined in Theorem [51) It follows that i?? = 1 — 
-Bi — i?2 = 1 — A2 is a projector and hence satisfies the optimality condition in 
Eq. (j25bp . Thus the measurement is optimal if and only if Eq. (|25ap holds, i.e., 

Ai£;?(72 - 7i)£;?Ai > 0. (30) 

Let us now assume that supp 71 n supp 72 ~ {0}. Then (1 — A2)Ai^ = 71^ 
and we arrive at the following 

Proposition 12. Let £ = {Ei, E2, E-?) be an optimal USD measurement for a 
pair of weighted density operators S = (71,72) with supp 71 n supp 72 = {0}. 
Then ti:{Eiji) = if and only 1/71(72 — 71)71 > 0. 

In this case the success probability is given by Psucci£',S) = tr(A272), and if 
£ is proper, then £ = (0, A2, 1 — A2). (A2 is the projector onto ker7i n supp 5. j 

Proof. Assume, that £ is optimal and satisfies tr(i?i7i) = 0. Then also for 
the corresponding proper measurement £' = {E[, E2, E!,) (cf. Sec. 12. 2p we have 
tr(iJj7i) = and hence 71(72 — 71)71 > follows. For the contrary, we have 
already shown that the if 71(72 — 71)71 > holds, then the proper measurement 
f = (0, A2, 1 — A2) is an optimal measurement. But due to Proposition fTT| this 
is the only optimal and proper measurement. Let now £ — {Ei, E2, E-r) be some 
optimal measurement, that is not proper. Then if the projection Ei of Ei onto 
supp 5 satisfies tr(i?i7i) = 0, then necessarily also tr(i?i7i) = holds. □ 
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Let us consider the situation, where the success probabihty for the states 
pi and p2 (both having unit trace) is analyzed in dependence of the a priori 
probabihty < pi < 1 of the state pi, while the a priori probability of p2 is 
P2 = 1 — pi. Then the optimality condition in Proposition [l^] is satisfied, if and 
only if for any G supppi, 

{l~ Pi){ip\p2\ip) >pi{Lp\pi\ip). (31) 

If there exists a \ip) £ supppi n kerp2 with \Lp) ^ 0, then this condition cannot 
be satisfied for any pi > 0. But if we assume supp pi nkerp2 — {0}, single state 
detection of P2P2 is optimal if and only if < pi < ^1, where £1 is given by 
(with \ip) G supppi and {ip\ip) = 1) 

^,^min/ , 1 = ^, (32) 

M \{^\ipi + P2Mj 1 + Ai' 

where {^/pi" denotes the inverse of on its support) 

Xi = mm{ip\^^ P2^^\f). (33) 

I'-p) 

The minimum in the expression for Ai is given by the smallest non-vanishing 
eigenvalue of the operator yfp\~ P2\fp\~ (remember, that we assumed supp p\ n 
ker p2 — {0}). Note that Ai > and hence there always exists a finite parameter 
range for pi, where single state detection of 72 is optimal. 

An analogous construction yields ^2, such that single state detection of 71 is 
optimal if and only if < p2 < ^2 • 



5.2 Fidelity form measurement 

An upper bound on the optimal success probability of USD was constructed 
by Rudolph, Spekkens and Turner in Ref. [5]. Let |7/i)(7/i| be a purification 
[201 m] of 7^, i.e., a positive semi-definite operator of rankl acting on an ex- 
tended Hilbert space M' ® J^huk, such that the partial trace over J^ux yields 
back the original weighted density operator, traux l7At)(7At| = In- Since the par- 
tial trace can be implemented by physical means, the optimal unambiguous 
discrimination of S = (71,72) cannot have a higher success probability than 
^pur _ (|-)/^)('^j|^ |72)(72|)- But is a pair of pure states, for which the op- 
timal success probability is known due to the result by Jaeger and Shimony 
[Hj. The map from S to 5^""^, on the other hand, can only be performed physi- 
cally in very special situations [321 US] ^^nd hence the success probability of 5^"'' 
in general only yields an upper bound. This bound is strongly related to the 
Uhlmann fidelity tr l^/pT^/pil of pi = 71/ tr(7i) and p2 = 72/ tr(72) [HI dS]. 
The Uhlmann fidelity is the largest overlap between any purification of both 
states pi and p2 ■ Due to this relation the bound was named fidelity bound |13j . 
In Ref. [ini dlj , necessary and sufficient conditions for the fidelity bound to be 
optimal where shown and the optimal measurement was constructed. In this 
section we summarize and extend these results. 

We continue to assume supp 71 H supp 72 ~ {0}. Herzog and Bergou showed 
in Ref. that the fidelity bound can be reached only if i??(72 — 7i)i?? = 0. 
(From Corollary m it is obvious that any such measurement is optimal.) Then 
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due to Eq. ^ wc find that any measurement with £^7(72 — 7i)i?v = is given 

by 

£■? = + (71 +72)"{7i72 + 7271 + \/7i-FiV7i + \/72-^2\/72}(7i +72)" 
= 1 - (71 + 72)"{%/7i(7i - -^"1)^71 + V72 (72 - ^2)V72}(7i + 72)", 

(34) 

where we abbreviated Fi = yj ^/ti 72-^71 and F2 = \J y/j2li \/l2- The converse 
is also true: 

Lemma 13. Let S ~ (71,72) be a pair of weighted density operators with 
supp7i n supp72 = {0} and let £ = {Ei, E2, E-?) be a proper USD measure- 
ment of S. Then i?772-E? — E-zjiE? if and only if E-? is given by Eq. (I34p . 

Proof. It remains to show the "if" part. First we multiply the identity 

(71 +72) (71 +72)^71 = 71 (35) 

from left by Qi as defined in the proof of Proposition [31 i.e., Qi is the oblique 
projector from ker72 fl suppiS to supp7i. We then obtain due to Lemma [TSl 
(cf. Appendix [X]) that 71(71 + 72)~7i = 7i and thus 72(71 + 72) 7i = 0. (One 
can show that 71(71 + 72)^ = Qi-) From the polar decomposition of -^/tTv^, 
it furthermore follows, that there exists a unitary transformation U, such that 
V7^^/72C/ = Fl (and hence [7^71 ^tJ = F2, cf. also Ref. [H]). Thus 

£'?V7i = (71 + 72)~{0 + 72\/7i + \/7i(\/7iV72^) + 0} 

= (71 +72)"{x/72(C/V7i/72)'^ +7iV72}C/ (36) 
= E7^U, 

i.e., we have E-?jiE-? = Ei^2E->. □ 

We refer to the measurement characterized by Lemma [T3] as fidelity form 
measurement due to the appearance of the operators Fi and F2, which satisfy 
l%/7i-\/72l = tri^i — trF2. According to the classification in Sec. 13.21 the 
fidelity form measurements are a (strict) superset of the measurement class 
[r, r], where r = rank 7172. (This can be seen from Lemma fTOl in the class [r, r] 
we have dimker(l — E-;) = and hence in particular -E?(72 — 71)^^? = 0.) 

Unfortunately, it is very rare that the operator given by Eq. is part 
of a valid USD measurement. The following Proposition states necessary and 
sufficient criteria. 

Proposition 14 (cf. Theorem 4 in Ref. p4j)- Let S = (71,72) be a pair of 
weighted density operators with supp7i fl supp72 = {0}. Then there exists a 
proper USD measurement £ — {Ei, E2, E?) of S with E-? given by Eq. p4p if 
and only if 71 - a/TtH^VTi > ^"'^ 72 - vV^TiVt? > 0. 

// the measurement exists, it is optimal and the success probability is given 
by Psucci£,S) = tr(7i + 72) - 2 tr ITtTVtJI- 

(Note that tr(7i + 72) — 2 tr | .^/ti .^72 1 is the square of the Bures distance 
[261 El [25] between 71 and 72.) 
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Proof. Due to the properties shown in the proof of Lemma 1131 for E-? given by 
Eq. we have 7i(l— £^7)72 = and ^/7^(1— £7)^/7^ = jf^—Ff^. Furthermore, 
one can write £'? = IIj^ + AA\ with (cf. Ref. [Ml) 

(71 +72)"(V7r+ V7?C/)/f\, (37) 

where U is a unitary transformation originating from the polar decomposition 
\fl\\fyiU = Fi (cf . proof of Lemma [13]) . 

Due to these properties, the necessary and sufficient conditions on E-/ shown 
in Proposition [3] reduce to the assertion of the current Proposition. □ 

If the criterion in Proposition 1 141 is not satisfied, then the optimal measure- 
ment cannot be of the form as given by Eq. (|34p . Thus by virtue of Lemma [T3l 
£'?(72 — 7i)-E'7 7^ holds and using Lemma ITOl we have that ker(l — E-/) 3 ker5. 
But due to Eq. it follows that ker(l — E'>) contains at least one vector 
either in supp7i or in supp72 (cf. Corollary 1 in Ref. [H]). 

Similar to the discussion of the single state detection measurement in Sec. 15. 11 
we ask for the values of the a priori probability < pi < 1 of pi , for which 
the fidelity form measurement is optimal. The first condition in Proposition [141 
71 — -Fi > 0, is satisfied if and only if for any \lp) £ supp7i. 



Pi{ip\pi\^) > ^pi{l-pi){ip\Ri\ip) (38) 



holds, where we abbreviated Ri ~ y/ ^/PiP2^/pi- Thus 71 — -Fi > if and only 
if TOi < pi < 1, where mi is given by 



mi = (39) 

with fii the maximal eigenvalue of y/pi ~ Ri^/pi ~ . With an analogous construc- 
tion we get that 72 — -F2 > if and only if m2 < P2 = 1 — Pi < 1- Then 

mi < pi < 1 - m2 (40) 

is the region where the fidelity form measurement is optimal. Note, that this 
region is empty when mi -I- m2 > 1. 

In summary, single state detection is optimal if and only if {(71 — F^ < 0) or 
(7I — F2 < 0)}, while the fidelity form measurement is optimal if and only 
if {(71 — Fi > 0) and (72 — F2 > 0}. The situations, where the optimal 
measurement is neither a single state detection measurement nor a fidelity form 
measurement seems to be related to the gap between "A > i?" and "A^ < 
for positive operators A and B. In the pure state case, however, A and B are 
of rank 1 and hence this gap does not exist. Indeed, in the pure state case, 
Ml = /^2 = tr(/3ip2) = Ai = A2 (where and £^ were defined at the end of 
Sec. l5.ip . Thus £1 = mi and £2 ~ ^2 and hence either the single state detection 
or the fidelity form measurement is always optimal. This is exactly the solution 
for the pure state case, as given by Jaeger and Shimony in Ref. [B]. 



6 Solution in four dimensions 

In this section we reduce the candidates for an optimal and proper USD measure- 
ment for the case where dimsuppiS = 4 to a finite number. These candidates 
are obtained by finding the real roots of a high-order polynomial. 
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Table 1: Possible ranks of the measurement operators of the optimal USD mea- 
surement in four dimensions (cf. also Sec. 13. 2p . 



rank Ei 


rank E2 


rank E-/ 


dimker(l - E-,) 


class 


cf. Sec. 
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[0,2] 
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2 
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2 





[2,2] 
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[1,1] 




2 
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2 


1 


[1,2] 




2 





2 


2 


[0,2] 
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Due to Proposition [2] and Proposition [51 it is sufficient to discuss the case of 
strictly skew pairs S — (71,72) with keriS — {0}, i.e., 

supp7i n SUPP72 = {0}, ker7i n ker72 = {0}, ^^^^ 
supp7i n ker72 — {0}, and ker7i n supp72 = {0}. 

Then, following the discussion in Sec. 13.21 there are six types of optimal USD 
measurements which differ in the rank of the measurement operators Ei and 
E2 ■ We list these possible types in Table [1] 

The measurements of the class [0, 2] and the class [2, 2] where already ex- 
tensively discussed in Sec. 15. H and Sec. 15.21 respectively. For the class [1, 2], the 
kernel of (1 — £'?) is one-dimensional. We will consider this class in Sec. 16.11 In 
the remaining case, where ranki?i = 1 = ranki?2, the measurement is a von- 
Neumann measurement (cf. Sec. I3.2|) . An example of this kind of measurement 
was first found in Ref. [TB]. Sec. 16.21 is devoted for a general treatment of this 
class. 

6.1 The measurement class [1,2] 

In this class the dimension of ker(l — £'?) is 1. According to Eq. IT^ . this kernel 
is either contained in supp7i or in supp72. Here we focus on ker(l — E-?) C 
supp7i, i.e., to the measurement type (1, 2); the case of ker(l — C supp72 
follows along same lines. Thus there exists an orthonormal basis of 
supp7i, such that 

E, = ,.\n){n\ + \(bM, (42) 
with \n) some normalized vector, orthogonal to ]0) and < < 1. 

6.1.1 The necessary and sufficient conditions 

Due to Proposition [21 the operator in Eq. (|^^ is a valid USD measurement, if 
and only if 71(1 — £'7)72 = 0, i.e., 

-fi\(t>^){(t>^\j2 = i^7i|"-)("-l72 (43) 

holds. This is equivalent to 

7iln) = a7il0^) and 72|rt) = &72|'/'^)- (44) 
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where a ~ {4>'^\n) and b — {h'{n\(p-^))^^ . Remember, that due to Theorem|3]we 
have \n) ^ ker7i, i.e., {n\(j)-^) ^ 0. On the other hand, from z/ < 1 it follows 
ah* > 1 and thus \n) ^ spand^-"-)}, i.e., \n) ^ supp7i. Given l^"""), we choose 
the phase of |n) such that {n\(j)^) > 0. 

The conditions in Corollary [9] can now be reduced to simple relations. By 
multiplying Eq. (|25bp from the left with |(/)-'-)(0-'-|(Ai - A2)~^, this first leads to 

(n|72-7i|n> =0. (45) 

Now it is straightforward to see that the conditions in Eq. (pSj) are equivalent 
to 

-^7(72 -7i)|n) = 0, and (46) 

(</>l72-7i|'/') >0. (47) 

By virtue of Eq. we can re-express Eq. in terms of |0) and 

yielding 

y^(<^^|7i|0^)(0^|72|0) = V^(</>^|72|0^)(0^l7i|'^), and (48) 

(49) 



a_ (</.^|72|</>^) 



The last equation enables us to construct z^|n)(n| from 

V^\n) = Vv{ji + 72)^^71 + l2)\n) 

= Vi^(7i +72)^^(071 +^72)10^) (50) 

where in the last step we used V abv = 1 and we abbreviated 

i^=(7i+72)"^(Va7&7i + yVa72)- (51) 

But due to Eq. \fajh is given in terms of |0) and Note, that if has 

full rank and hence ensures v = {(t)^\K^ K\(f)^) > 0. 

We summarize: The optimal USD measurement is of type (1, 2) if and only if 
there exists an orthonormal basis (10), 14'^)) o/supp7i, such that the conditions 
in Eq. (gT]), in Eq. (gS]), and {cfi^lK'^ K\(t)^) < 1 are satisfied. 



6.1.2 Construction of a finite number of candidates for E? 

In the following we will show, that already Eq. ((48|) reduces the possible can- 
didates of span{|0)} to a finite number and hence the remaining positivity 
conditions can be easily checked. Eq. is a complex equation. Thus the 
absolute value and the phase of the left hand side and the right hand side have 
to be identical. This leads to 

(0|72|0^)(0^|72|0)(0^|7i|0^> = ('^I7i|'^^>(0^l7i|0>('^^l72|0^), (52) 
{4^hM{cl)\li\cl)^)>0. (53) 
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Let (|si), |s2)) be an orthonormal basis of supp7i, such that 7i|si) — gii\si) 
with > .912 > 0. We abbreviate g2i = (si|72|si) > 0, = g^i - g^2, and 
523 — (si|72|s2)- We ensure 323 > by choosing a proper global phase of |s2)- 
In the case gi — we use a basis where 1723 = 0. 

First we consider, whether \si) or |s2) is a candidate for |(/)). In either case, 
Eq. can only be satisfied, if 523 = 0. From Eq. (|T7)) we find that \si) is a 
candidate only if 321 > ffii and analogously, \s2) only if 522 > .912- 

We now assume that neither of the above two cases is optimal. Then any of 
the remaining bases (apart from global phases) are parametrized by 

10) = (l + x-2)-^(|si)+xe^''|s2)), (54a) 
\cl>^) = {l + x^)-Hxe-^'>\s,)-\s2)), (54b) 

where x 7^ is real and —§<"!?<§■ 

Using these definitions, Eq. can be written as 

ffi523 sini9 = and xgi{xg2 + 523(2^^ - 1)) > 0. (55) 

Let us first discuss the case where 523 7^ 0. In this situation we geiH from 
Eq. (|55p the phase = 0. The solutions of Eq. ([5^ are given by the real roots 
of the polynomial of degree six in x, 

x^glix^g^i + .922 - 2x523) - {xg2 + [x^ - l)g2if{x^gii + 512) = 0. (56) 

(Since 5I3512 > 0, this polynomial cannot be trivial.) 

It now remains to consider the special case, where 1723 = 0, but neither 
10) — \si) nor 10) — \s2) is optimal. If 523 = 0, then Eq. (l52|) reads 

x^i.glg2i - glgii) = 5I512 - g1g22- (57) 

Assume that this equation has some solutions where x is real (there might be 
infinitely many). None of these solutions leads to an optimal measurement, 
as we exclude by the following argumentation. Neither Eq. (|T7|) nor Eq. 
depend on From the necessary and sufficient conditions for optimality (cf. 
end of Sec. [gXT|) . only the condition ly = {(l>^\K^ K\(l>^) < 1 remains to be 
satisfied. Since Eq. (|49p does not depend on d, also K is independent of -d. 
Thus v as a function of ?? is of the form 

= (1 + a;2)-i (x^ku - 2 xReie-'"^ k^) + K22) , (58) 

where we defined K.y = {si\K^ K\sj). In particular this function is continuous 
in 'd. Assume, for a given x 7^ (satisfying Eq. (I47|) and Eq. (|48|) ). there exists 
some —tt/2 < < tt/2, such that i'{'d) < 1. Then there exists an e > 0, such 
that also i) + e < Tr/2 and i'{-d + e) < 1. Hence for i) as well as i? + e all optimality 
conditions would be satisfied. But for both values, we get a different vector |0) 
and hence there would be two different operators E-?, both being optimal. This 
is a contradiction to Proposition [TTl It follows that for 923 = 0, only |0) = \si) 
and 10) = \s2) can yield an optimal solution. 

^Remember, that we chose our basis such that 323 ^ implies gi 7^ 0. 
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6.1.3 Summary for measurement type (1,2) 

Let us briefly summarize the algorithm to obtain the optimal measurement E-? 
for the case where ranki^i — 1 and rank £'2 = 2. 

We construct some basis (|si), |s2)) of supp7i as described in the paragraph 
below Eq. ([55)1 . In a next step, we construct candidates for the basis (|(/)), \4>'^))- 
There are two cases: 

(i) .923 = 0. If 321 > .911, then = |si),|0-'^) = |s2)) is a candidate. If 
522 > 512, then (10) — |s2), 10^) — \si)) is a candidate. 

(ii) 523 7^ 0. For any real root a; 7^ of the polynomial in Eq. ([56|l . the basis 
(10), 10-L)) as defined in Eq. ^ is a candidate, where -d = 0. A candidate 
in addition has to satisfy the second part of Eq. (|55p and Eq. (|T7)) . 

For any of the candidate bases (if any), we construct v^l") using Eq. ([50|) . 
At most one of the bases will satisfy ly = Y^(n|n)y^ < 1. If such a basis exists, 
the optimal measurement is given by Eq. (|42p . 

6.2 The measurement class [1,1] 

In Sec. 13.21 we have already seen, that if rank £^1 + rank i?2 = rank7i(= rank 72), 
then both Ei and E2 are projectors. Hence there are orthonormal bases (|'0i) : li'i)) 
of ker72 and (|'i/'2), IV'2")) of ker7i such that 

^i^lV-iX^il and £;2 = IM^2|. (59) 

Since 1 — £1 — £2 > 0, necessarily (V'ilV'2) = must hold. Using this notation, 
Eq. (|^ is equivalent to 

\{^^\^2)f{HlM > (V'^bilV'^) (60a) 

|(V'2'-|V'i)l'(V'i|7i|V'i> > (^2'-|72|^2^) (60b) 

(^2^l^i)(^i|7i|^i^> = (^2 IV-^) (^2^ 1721^2), (60c) 

while Eq. (j25bp is satisfied identically. (Note that these equations only follow if 
all vectors are normalized.) 

6.2.1 Construction of a finite number of candidates for E7 

Let {\kii)) and {\k2i)) be Jordan bases (cf. Appendix [C)) of ker7i and ker72, i.e., 
{\kii)) and (|fc2i)) are orthonormal bases of ker7i and ker72, respectively, such 
that (fci/|A;2;) > and for i ^ j one has {kii\k2j) = 0. Due to our assumptions in 
Eq. (|4T|) we have < {kii\k2i) < 1. We choose (A:ii|fc2i) > (fci2|fe2)- In case of 
degenerate Jordan angles (i.e., (A:ii|fc2i) = (^121^22)), these bases are not unique 
and we then choose bases, such that (A:2i|7i|A:22) = 0. We abbreviate 

gii = {k2ihi\k2i), g2i = {kii\j2\kii) , 

513 = I(fel|7l|fc22)|, 523 = |(fcll|72|fcl2)|, (61) 

and choose the global phases of the vectors \ki2) and |fc22) in such a way that 
still (/ci2|fc22) > but also 5136*"^ (A:2i|7i|^22) and 5236"**^ = (fcii|72|fci2). 
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with < ip < TT. We let = if (713 = or 323 = 0. Finally we define 
c = {ki2\k22) / {kii\k2i) ■ 

After choosing the Jordan bases, let us first consider the case where = 
\k2i). Then one has (fixing the phases) |i/'2) — \ki2), \ipi) — 1^22), and IV'2") = 
|fcii). Eq. (I60cp reduces to C523 = 513 and sinip = 0. Analogously, in the case 
iV'i) = 1^2), we obtain cgi3 = (723 and sin ip = 0. 

We now consider the case, where is not one of the basis vectors |fc2j)- 
We fix the global phases of {ip^) and and choose the parametrization 

IV'i) = (l + -^')"^(|fc2i)+:cc"'|fc22)), 
IV'^> = (l + -^')-^(^e-''|fc2i)-|fc22», 

1^2) = (1 + c^xYH-xe'^'^clkn) + |fci2)), 
1^2^) = (l + c2x2)-^(-|fcn) -xe*''c|fci2)), 

with the real parameters x ^ and —tt/2 < d < tt/2. 
Using these definitions, Eq. (|60cp reads 

ix'c' + If [xg^ - e^(''+^)5i3 + x'e-^^^+^^ g,,) = 

c{x' + 1)2 [xcg2 - e^(''-^).g23 + x^c'e^^^"-^^ 923) (63) 
The imaginary part of this equation gives 

{x^c^ + l){l + x^) 

■ [{x^c' + 1) sin{d + (p)5i3 - (1 + x^) sm{d - ^)cg23] = 0, (64) 

which can only hold if already the term in square brackets is zero. This is the 
case if and only if 

sini? = A2 COS!?, (65) 

where 

Ai = (c(a;2 + 1)523 - {x'^c^ + 1)313) cos(p, (66a) 

A2 = (c(a;2 + 1)523 + {x'^c^ + 1)313) sin^. (66b) 

In order to get the solutions of Eq. (|63p . we consider now its real part, 

{x^c^ + l)^{xgi + cos(i? + (p)gi3{x^ - 1)) = 

c(x2 + l)2(a;c52 + cos(i? - ip)g23{x^c^ - !))■ (67) 

Using the abbreviations 

Bi = (x^c^ + 1)^x51 - c{x'^ + \ fxcg2, 

B2 = [{x'c^ + lfgi3{x'' - 1) - c{x^ + lfg23{x^c^ - 1)] cos 9?, (68) 
B3 = [{x^c^ + l)'ffi3(a;' - 1) + c(x2 + l fg23ix^c^ - 1)] sin^, 

we get the equivalent expression 

Bi ^ B3sm'd ~ B2Cos-d. (69) 
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Taking the square of this equation, multiplied by {Af + A^), we obtain due to 
Eq. 1^5]) the polynomial (with a degree of at most 8 in a;^) 



BliAl + Al) - {A1B2 - A2B^f = 0. (70) 

This polynomial is trivial if and only if (713 — 523 = 0. (In order to see 
this, we consider the highest and lowest order term, which only can vanish, 
if -(cgi3 - 523)^cos^(/9 = (cgi3 + g2zYsiv? ip and if -(0^23 - 5i3)^cos^(p = 
(c(723 + <?13)^ sin^ Lp, respectively. But due to our particular choice of the bases, 
this can only hold if already (713 = (723 = 0.) It is straightforward to see, that in 
this case none of the conditions in Eq. depend on Now suppose there is 
a solution of these conditions with a; ^ 0. Then any possible value of d leads to 
a different, but optimal measurement, in contradiction to Proposition 1111 

In any other situation we get from Eq. (|70p a finite number of real solutions 
X ^ Q. The corresponding value for can be obtained as follows. If Ai ^ 0, 
then from Eq. fSK]) we have -d = arctan(A2/Ai), while if ^1 = and A2 7^ 0, 
then d = — 7r/2. If Ai = A2 ~ 0, then soup = and 

^2^l£52^^^ (71) 
c cgiz - 323 

where from > it follows that {cg2z — 5i3)(c5i3 — 323) > 0. From Eq. ((69l) 
we have 

2 cos 1? 313323 (cg23 - 513) = xc{gl^gi - 51332), (72) 
which can be used in order to obtain i9. 



6.2.2 Summary for measurement class [1, 1] 

In a first step we construct Jordan bases as described in the first paragraph of 
Sec. 16.2.1] Then we collect candidates for the bases 

(i) If sin(^ = and C523 = 5i3, then (IV^i) = jfei), IV"]^) = [^22)) and (IV'2) = 
1^12), 1^2^) = 1^11)) is a candidate. 

(ii) If sin(^ = and C313 = 323, then (j-f/'i) = 1^2), \'4'i) = |fei)) and (IV'2) = 
1^11), = 1^12)) is a candidate. 

(iii) For any real root x 7^ of Eq. (|70p. we get a unique value for d from 
Eq. jM]) (if Al 7^ or ^2 7^ 0) or from Eq. jH]) (if Ai = and A2 = 0). If 
Al 7^ or ^2 7^ 0, then in addition Eq. ([69|) has to hold. For each [x, t?), 
we obtain the candidates from Eq. (|62p. 

At most one of the candidates will satisfy Eq. (j60a|) and Eq. (j60b|) . If such 
a candidate exists, the optimal measurement is provided by Eq. (|59p . 



6.3 Examples 

We want to discuss a few examples of USD in four dimensions, which demon- 
strate the structure of the previous results. We consider three examples which 
belong to case [vii] in the flowchart in Fig. [31 Thus these pairs of states in par- 
ticular are strictly skew and do not posses a common block diagonal structure. 
In Fig. [T] the optimal success probability of the two states 
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Figure 1: Optimal success probability of the states given in Eq. ([75)) . depending 
on the relative probability pi of the occurrence of pi (solid line). The dashed 
lines denote the success probability of single state detection (lower bound) and 
the dotted line corresponds to a simple upper bound. In brackets we denote the 
measurement types as defined in Sec. 13.21 
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is given in dependence of the a priori probability of pi of pi (solid line). Follow- 
ing the results of Sec. 15.11 and Sec. 15.21 we can directly calculate the probability 
range where single state detection (class [0,2]) and the fidelity form measure- 
ment are optimal (roughly class [2,2]). The remaining optimal measurements 
belong to the classes [1,2] or [1,1] which are obtained by following Sec. 16.11 
or Sec. 16. 2[ respectively. We also plotted the "bound triangle" which can be 
easily calculated for any pair of density operators. The lower bounds corre- 
spond to single state detection (dashed lines) and the upper bound connects 
the points where single state detection stops being optimal (dotted line). The 
latter one is an upper bound due to the convexity of the success probability 
function Psucc(pi). 

The next example (Fig. [2] (A)) was found by generating "randomly" pairs of 
density operators. The data of the states is available as supplemental material. 
This example shows that the measurement types (1,1), (1,2), and (2,1) can 
appear in more than one probability range. This is not possible for the other 
measurement types. It is also an example for a pair of states, where all possible 
measurement types can be optimal in some probability range. 
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Figure 2: Optimal success probability depending on the relative probability 
Pi of the occurrence of p'l- (A): A "random" example. (B): An asymmetric 
example, cf. Eq. (|7^ . In brackets we denote the measurement types as defined 
in Sec. 13.21 For details see text in Sec. 16.31 
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Yes J Pure state USD 

\ in each block.|7^ 



Measurement 
class [0,r] or [r,r] 



Measurement 
class [1,2] or [1,1] 




Check optimality 
of bounds, [vili 



Figure 3: Flow chart of a generic strategy to solve/reduce the USD of two mixed 
states. For details see text in Sec. [7] 

The last example (Fig. [D (B)), given by 



Pi 



5 41 

P'2 = + ^I^X"^! 



im 



with 



\w) = ^ ((-1)1/7 ((722 + 275)10) + (V22 - 2^/5)|l)) + 2VTQ (|2) + |3))) , 
\v) = (((-1)^/^1)10) + + 2V2(-1)1/5|2)) , 

(74) 

is devoted to show that there is no deeper general structure which optimal mea- 
surement classes can be connected with each other. Here, the fidelity form mea- 
surement directly follows the single state detection measurement. This example 
also shows a high asymmetry in the success probability function Psucdpi)- 



7 General strategy 

In this section we want to summarize the known results for the unambiguous 
discrimination of two mixed states by suggesting a strategy in order to find the 
optimal success probability, cf. Fig. [31 

In step [i] we check whether the pair of weighted density operators S — 
(71 7 72) is strictly skew. A pair S is called strictly skew, if supp7i n supp72 — 
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{0}, supp7i n ker72 = {0}, and ker7i n supp72 — {0}. If the pair fails to 
be strictly skew, one can use the reduction method described in the discussion 
of Lemma [7] (Fig. [3] [ii] ) . After this reduction, it suffices to find the optimal 
measurement £,' of the strictly skew pair S' . 

From now on we assume that the states are strictly skew. In the next step 
(Fig.[3][iii]) we check with the commutator criteria shown in Ref. |19j whether the 
two states have an at most two-dimensional common block diagonal structure. 
If this is the case in Ref. [I^ it was shown how to construct diagonalizing Jordan 
bases, which then can be used to compose the optimal measurement from the 
pure state case (Fig. |3] [iv]), as shown in Ref. [l8lfT9]. 

For states without an at most two-dimensional common block diagonal struc- 
ture, optimality of one of the two generally solved measurement classes, i.e., 
single state detection (Proposition[T21) or the fidelity form measurement (Propo- 
sition [T3|), can be checked (Fig. [3] [v]). Otherwise the optimal measurement can 
be calculated if the two states effectively act on a four-dimensional Hilbert-space 
(dim suppiS = 4). Here the optimal measurement of the two remaining measure- 
ment classes [1, 2] and [1, 1] can be calculated and checked according to Sec. 16.11 
and Sec. 16.21 respectively (Fig. [3] [vi]). In principle one could also find optimal 
measurements of states which have a four-dimensional common block diagonal 
structure, because then blockwise combination of optimal measurements yields 
the optimal measurement for the complete states. Unfortunately, there is so far 
no known constructive method to identify such blocks. This question is left for 
further investigations. 

The last possibility is to check optimality of upper and lower bounds on 
the optimal success probability (Fig. [3] [viii]). Examples of such bounds were 
presented e.g. in Ref. [H |T31 [HI [27]. For some of these bounds, optimahty 
conditions are known, while in the remaining cases one can use Corollary [9] in 
order to check for optimality. 

If the procedure sketched here fails to deliver the optimal solution then there 
is up to now no systematic method known to find an analytic expression for the 
optimal USD measurement. 

8 Conclusions 

We analyzed the unambiguous discrimination of two mixed states p\ and p2 and 
the properties of an optimal measurement strategy. We first showed (cf. Propo- 
sition [3]) that any unambiguous state discrimination measurement is completely 
determined by the measurement operator E-i^ the operator which corresponds 
to the inconclusive measurement result. A further analysis for optimal mea- 
surements showed (cf. Theorem Hj) that the rank of E-i is determined by the 
structure of the input states, ranki?? = rankpip2- 

This fact leads to one of our main results, namely, that the optimal measure- 
ment for a given pair of mixed states and given a 'priori probabilities is unique 
(cf. Proposition [11]) . This uniqueness might have interesting consequences, e.g. 
in the analysis of the "complexity" of the optimal measurement. Interest- 
ing questions here are whether the optimal measurement is von-Neumann or 
whether the optimal measurement is non-local, as e.g. discussed in Ref. |28j . 

Eldar et al. in Ref. jl7j provided necessary and sufficient conditions for a 
measurement to be optimal, but in many situations this result was not opera- 



26 



tional. Wc simplified these optimality conditions in Corollary[51 which now can 
be appHed directly to a given measurement. 

As an application of this result, we analyzed the single state detection case, 
where the measurement only detects one of the two states. Although this mea- 
surement may seem to be pathological, it turns out to be always optimal for a 
finite region of the a priori probabilities of the states. We derive an analytical 
expression for the bounds of this region. 

Finally, we constructed the optimal measurement for the unambiguous dis- 
crimination of two mixed states having rank pi = 2 = rank p2 (due to a results 
by Raynal et al. in Ref. [11] , this construction can be extended to the case where 
one of the density operators has rank 2 and the other state has arbitrary rank) . 
The solution splits into 6 different types, each of them requiring a different 
treatment. This solution is analytical, but in certain cases the roots of a high 
order (up to degree 8) polynomial are needed. Due to the complicated structure 
of this solution it may turn out to be quite difficult to analyze the next step, 
namely the discrimination of two mixed states with rankpi = 3 = rank/92. 

It would be interesting to find strategies in order to detect symmetries and 
solvable substructures in unambiguous state discrimination, e.g. to find four- 
dimensional common block diagonal structures, analogously to the result for 
two-dimensional common block diagonal structures in Ref. [19\ Also, there are 
only a few results on the optimal unambiguous discrimination on more than two 
states [221 [301 1311 1321 [IZl 133] • We think that several concepts presented in this 
contribution may also generalize to the discrimination of many states. 
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A Appendix: Technical statements 

We put the following two Lemmata for completeness, but without a proof. 

Lemma 15. Let jz/, ^ and he subspaces of M' with jz/ _L and SS . 
Then si ^{SS + = si ^ SS . 

Lemma 16. Let A and B he operators on with BJif ker A — {0}. Then 
kei AB = kerS. 

The useful Eq. (fT2|) is based on the following 

Lemma 17. Let X he an operator and (11^) he a family of projectors with 
ker(Efenfc) = {0}. For k^l, assume that HfeXn, = 0. (Note that UkU-i ^ 
is allowed.). Then ker A" = '^k^^''^ ^ ^k-J^)- 

Proof. The "d" part is obvious. For the contrary let |$) e ker X. Then there 
exist vectors \Lpk) G Hfe^ such that |$) = J^k Iv'fc)- We have 

J2T^iX\iPk) ^IikX\ipk) ^llkX\<^^) ^0. (75) 
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Since J2k^k has a trivial kernel, it follows that X\(fk) — 0, i.e., \(pk) G kerX. 

□ 

The following Lemma lists properties of non-orthogonal projectors, i.e., idem- 
potent operators which are not self-adjoint. Such operators where considered 
e.g. in Ref. [SH El]. Each of the statements (ii) - (iv) is a valid definition 
for the oblique projector, where (ii) is the formal definition, (iii) is an explicit 
construction and (iv) are the central properties for our purposes. 

Lemma 18. Let A and II be two (self- adjoint) projectors on Jif with AJif D 
(IlJif)-^ = {0} and (A^)-^ n = {0}. For an operator Q the following 
statements are equivalent: 

(i) Q is the (unique) oblique projector from Ajff to IlJif. 

(ii) g2 = Q, Q.ye = njr, and kerg = kcr a. 

(iii) Q is the Moore-Penrose inverse of AH. 

(iv) QA ^Q,UQ ^Q, AQ = A and QU = U. 

Proof, (ii) formally defines (i). In order to show that (iii) follows from (ii), 
we write Q in its singular value decomposition, Q = J2 with qi > 
and (Itt^)) and (|Ai)) an appropriate pair of orthonormal and complete bases of 
nj^f and Ajf , respectively. Then Q'^ = Q is equivalent to {\k\nk) = q^^ and 
(AilTTj) = for j ^ j. Since AH = q~^\\i){TTi\, it immediately follows that Q 
is the Moore-Penrose inverse of AH. 

From the explicit form ofthe Moore-Penrose inverse, Q = J2i (Ai|7rj)^^|7rj)(Ai|, 
one easily verifies the properties in (iv). 

We have from (iv) that QJif = UQ.Jf C HJT = QUJf C QJf and hence 
Q.Jf = njf and analogously kerQ = (Q^'J^)^ = {AJ^)^ ^ kerA. Finally, 

Q2 = g(ng) = (gn)g = ng = g. □ 



B Appendix: Proof of Corollary [9] 

B.l Necessity 

Let us abbreviate p. ^ 3 — fi, ~ A^ZA^ and = A^ZA^. From Theorem [S] 
it follows that for any optimal measurement we have {fi — 1, 2) 

Z,, - A,,7^A^ > 0, (76a) 

(Z^ - A^7^A^)£;^ = 0, (76b) 

and 

Z^{A^^Ef,)^YpE-^A^, (77a) 

Y^,{Af,-E-^)^ Z^Ef^Af,, (77b) 

where the last two equations follow from A^ZE-rA^, = with E-r = 1 — E^ — Ep_. 
From Z > we have Z„ > and Yr, = Y.l We find 



Z^ - A^7^A^ = (A,, - E^)iZ,, - 7^)(A,, - E^) 

- A,,Ef,Y^iA^ - E,,) - {A,, - EMAf^ - 
^ A,,E-^-1-^Ef,A^ - {A,, ~ E^)-i^{A^ - E^) 
^ A^,E-ap,E-?A^ - A^E7j^E7A^. 



(78) 
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Thus Eq. (|^^ . Eq. (p6b)) follow from Eq. ((7^ . From Eq. ((76b)) we have 

A^£;?(7p-7p)^^?^^M = 0- (79) 
Combining Eq. (|77ap and Eq. (|77b|) we obtain 

FpA^ = Z^B^A^A^ + Z^(A^ - E^) (80) 
and hence from Eq. (I77a|) and due to Eq. (I76bp , 

ApZ^(A,, - Ef,) - (i;A^)t£;^A,, 

= Kp,Kf,E-^-i-^Ep,k^ + {Kf, - E^)jf,E^,Af, (81) 
- Ap(A^ - £;p)Z^(A^ - Ef,) + (Ap - £;^)7^£;^A^, 

where in the last step we used the result from Eq. (|78p . Thus we have found 
A^,Ef,jf,{Af, - £;^) = (A,-, - Ef,)jf,EfiA^,, i.e., Eq. follows. 

This equation together with Eq. ([7^ for /i = 1 and fi ~ 2, finally yields 
Eq. ([25b| . 

B.2 Sufficiency 

We first get rid of the non-skew parts of S: 

Lemma 19. With the definitions and preliminaries as in Proposition^^ if E-? 
is a proper USD measurement of S and satisfies Eq. (|25p for S, then Ef^°^ — 
i?7 + (1 — Ilskew) is a proper USD measurement of S^^™ and satisfies Eq. (^5]) 
for 3^^"=^. 

Hence, if we further show, that from the second part of the Lemma, it follows 
that Ef^'^'" is optimal for S^^'^'", then we have due to Proposition [HI that also 
E^? is optimal for S. 

Proof of Lemma\19[ We denote by Ei the projector onto supp7i n ker72, by 
T,2 the projector onto supp72 D ker7i and by 11^ the projector onto supp7^. 
Then multiplication of Eq. (|25ap by Ei yields 

i:iE'n2E?^i - J:iE-niE->T,i > 0. (82) 

Since for a USD measurement IliE->"/2 — ^ij2, we have S172 = and only the 
second term remains, which henceforth must vanish. This yields ^^Ef'Ei = 
or equivalently IIii^vEi = 0. A further multiplication from the left by Ei to- 
gether with the property E-? > proves that E'vEi ~ 0. An analogous argument 
can be used in order to show i??E2 = 0. Now, following the same lines of ar- 
gument as in the proof of Proposition [51 it follows that Ef^""" is a proper USD 
measurement for 5*^*^°*. 

From -EvE^ = it in particular follows that i??(72— 7i)£'? = £'f'^™nskcw(72 — 
7i)n,kcwi^?''™. Using 

ker nskow7pnskcw = ker7p(l - E^) = E^Jf © ker7^, (83) 

it is now straightforward to show that E^^™ satisfies Eq. ([^ for S^^™. □ 
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It remains to consider the case where supp 7inkcr 72 = {0} = supp 72nker 71 . 
We define i?^ to be the obhque projector from kei^p, to ker7^. Note that = 
Rp^. We furthermore denote by the obhque projector from ker7^ n supp 5 
to supp 7^ n (ker7i + ker72). Then the multiphcation of Eq. (|25bp by from 
the left and by from the right yields 

ApE-!{jp~j^)E,Ef, = 0. (84) 

Let us define 

Vf, ^ A^£;7(7p - ii,)E,Kp + A,,7^A^, (85) 

W-p = {Rf, ( ~ E-p) + KpE-p)Vj,. (86) 

Then, using Eq. (PBc)) . ((5I| . we have 

Vp{K^, - Ef,) = k^,E7{-if, - lt,)E-?Af, + Af,-/pE-?Af, 
= A^,E-?-/f,E-?Ap + Ep-/pE7Af, 
= -ApEf,jf,E7Ap - Rf,E->-ipE-?A^, 
= -{A^Ef, + R-pE7)j-pE7Ap 
- (ApEf, + i?p(A^ - E-p))-f-pEf,Ap 
= Wf,E-pA^. 

And a similar equation holds for Wp^: 

Wp{Ap - Ep) =RpE,hp - jp)E'!Ap + Rp{Ap ~ Ep)jp{Ap 
+ ApEpjpiAp ^ Ep) 
^RpE7jpE7Ap - ApE-npE7Ap 
= - i-Ep - (A^ - ii;^))7^£;^A^ 
^VpEpAp. 

We combine these two equations and find Wp — Vp{EpAp + (A^ — Ep)Rfj), i.e., 
in comparison with the definition in Eq. (|86|1 . Wp = Wj^. Furthermore, we have 
WiV-^Vi = Wi and V2 = WiVj-W}, where 1/f denotes the inverse of Vi on 
its support. (The second identity follows from WiVfW} = Wi(Ai — i?i)i?2 + 
W1E1A2.) 

We construct the operator Z as Z = TViT^ with T = Qi + Q2WiVf . Since 
Vi > we have Z > 0. From our previous considerations, A^ZA^ = and 
A^ZA^ = Wp follows. Then using Eq. ((26a|) . Eq. (l26bl) . and Eq. §^ it is 
now straightforward to verify that Z satisfies the conditions in Eq. (j24bp of 
Theorem [Sj 

It remains to show that ZE-? = 0. First, with Hy the projector onto supp 71 n 
supp72, we have ZE-tli]^]^ = ZHy = due to IInQ^ = 0. Analogously with IIj^ 
the projector onto ker5, ZE7ll± = ZII^ = 0. Thus we only need to show that 
ApZE7A^ = 0. But this follows from Eq. §7^) and Eq. jMI). 

C Appendix: Construction of Jordan bases 

In this Appendix an explicit construction of Jordan bases (also sometimes called 
canonical bases) of two subspaces and ^ of is given (cf. also Ref. [551 M 
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[inillH])- Jordan bases are orthonormal bases {\ai)) of £/ and {\bj)) of such 
that {ai\bj) — for all i ^ j and {ak\bk) > for all k. With 5*^ we denote 

the d X dimensional matrix where the columns are given by the basis 
vectors of some orthonormal basis of £/. Analogously we define and by 
using ^. 

Consider a singular value decomposition of 

SlSgg^U^DUl, i.e., iS^U^)HSsgU,j) ^ D, (89) 

where D is a Un/ x dimensional diagonal matrix. U,s/ and are unitary 
matrices. Let us denote the i-th column of SsiV by \ai) and the j-th column 
of SagUgg by \hj). Then (|ai)) and are Jordan bases of ^ and 
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